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o : 

' Abstract. We consider a Magnetic Rotator Wind-Disk(MRWD) model for 

I the formation of Keplerian disks around Be stars. Material from low latitudes 

i- of the stellar surface flows along magnetic flux tubes and passes through a shock 

O ' surface near the equatorial plane to form a pre-Keplerian disk region. Initially, 

m I the density in this region is small and the magnetic field helps to maintain 

. super-Keplerian rotation speeds. After a fill-up time, the density of the disk 

' is significantly larger and the magnetic force in this region becomes negligible 
compared with the centrifugal force. The material then moves outwards to form a 

T-H ■ quasi-steady Keplerian disk. During the fill-up stage, the meridional component 

^ ' Bm,* of the magnetic field at the stellar surface must be larger than a minimum 

OO , value Bm.*,min- Thc radial extent of the quasi-steady Keplerian region will be 

CO i larger when Bm ^, is larger or when viscosity plays a role. In B-type stars, the 

values of Bm * mm are of order 1 G to 10 G. We find that a condition for the 

, formation of shock-compressed disk regions is that in faster rotating stars, the 

■ wind speed must be correspondingly larger. 

o 

Pi 1. Disk Formation Processes 

. The MRWD (Magnetic Rotator Wind-Disk) model for the formation of Keple- 

^ ' rian disks involves two processes. The first is a fill-up process in which material 

•th , from a star flows along magnetic field lines and passes through a shock sur- 

[ face to form a pre-Keplerian disk region. The magnetic field supplies angular 

momentum to the wind and assists the flow of material towards the disk re- 
gion. Initially, the density in the pre-Keplerian region is small and it increases 
significantly during this process. 

The second process occurs after the magnetic force in the disk region be- 
comes small compared to the centrifugal force. Then, super-Keplerian material 
in the pre-Keplerian region expands to form a quasi-steady Keplerian disk. If 
there is no viscosity or inflow into the disk region, angular momentum will be 
conserved. The time at which the second process starts depends on the rota- 
tional speed, magnetic field strength and rate of fill-up of the disk region. In some 
stars, sub-Keplerian wind material may continue to flow into the super-Keplerian 
disk region. Also, as the angular velocity of the expanding disk approaches a 
Keplerian distribution, magnetorotational instability (MRI) and viscosity can 
become important. When there is inflow or viscosity, the radial extent of the 
quasi-steady Keplerian disk region will be larger. In situations where the mag- 
netic force in the disk is stronger than the centrifugal force at the end of the 
first process, the disk will continue to maintain its pre-Keplerian structure. 
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Figure ^ gives a schematic meridional cross-section of a pre-Keplerian disk 
during the fih-up stage. It shows the magnetic field lines from the stellar sur- 
face to the inner and outer end points Xinn and Xn^, respectively, of the pre- 
Keplerian region. In a meridional plane, streamlines of flow from the star to the 
disk coincide with magnetic field lines. The limiting streamline from Pum to Xum 
separates streamlines that flow into the disk region from the open streamlines 
that flow away from the star and the disk. Let X be a point on the equatorial 
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Figure 1 . Schematic meridional cross-section of a pre-Keplerian disk during 
the fill-up stage. 



plane in the pre-Keplerian region. Let P be the point on the stellar surface that 
is connected to X by a magnetic field line that loops back to a point P' on the 
star. If to is the angular velocity at P and P' , t he angular velocitv at X during 
the pre-Keplerian stage will also be equal to oj (|Maheswaranl 1200^ 1 . The ratio 
of the angular velocity at P to the critical angular velocity at the same point 
IS Q = a; 

i?3/2/ [GM{1 - r)]^/^ where M and R are the stellar mass and radius, 
respectively. T is the ratio of the continuum radiation force to gravity. For 
points on the equatorial plane we put x = r/R where (r, 6, (p) are spherical polar 
coordinates. Let X^ep and Xesc, respectively, be the points in the pre-Keplerian 
region at which the rotational velocities are equal to the Keplerian and escape 
velocities. Then, Xkep = a"^/^ and Xesc = 2^/^a~^/^. 

When the meridional component -Bm,* of the magnetic field at the point 
Piim on the stellar surface is increased, the value of xnm also increases. The 
MRWD model requires that a super-Keplerian region be present in the pre- 
Keplerian disk. Hence, Xum must be further away from the star than X^ep so 
that xiim > Xkep- Thus, Bm,-), must be larger than a minimum value Bm,*,min, 
which is the meridional field strength when Xum coincides with X^ep. Figure [21 
schematically shows the diff'erent possible locations of X^nn and X^^ in relation 
to Xkep and X^sc- During the second process, material between X^ep and X^sc 
will move into Keplerian orbits. When there is no viscosity or inflow of wind 
material, the material beyond X^gc will flow away. If viscosity or inflow becomes 
important, material in the region between X^sc and Xnm can become part of 
a quasi-steady Keplerian disk. A pre-Keplerian region having configuration 1 
in Figure [21 will yield a quasi-steady Keplerian disk region with maximal radial 
extent. The optimal model corresponds to configuration 5, where X^^ coincides 
with Xesc and Xinn coincides with X^ep or is closer to the stellar surface than 
Xkep- This model requires the least magnetic field strength at the stellar surface 
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to form a Keplerian disk with the largest radial extent in the absence of inflow 
or viscosity. 
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Figure 2. A schematic picture of the different possible locations of the 
end points Xinn, of pre-Keplerian disks in relation to the critical points 

The Keplerian region formed when i?m,* is only slightly larger than Bm,*,min 
will be in the shape of a ring. For larger values of B^,* the radial extent of 
the Keplerian region will be larger. In an inviscid model without inflow, the 
Keplerian region has its largest radial extent when -Bm,* equals Bm,-k,opt, which is 
the field strength required for the optimal model. For values of -Bm,* larger than 
Bm,*,opt there is no increase in the radial extent of the quasi-steady Keplerian 
region as long as the material beyond Xggc moves away from the disk region. 

When the rotation rate a of the sector PumPinn is approximately constant, 
the inner radius of a quasi-steady Keplerian disk is a~'^^^R. Let ^end be a point 
on the outer boundary of the quasi-steady Keplerian disk. When there is no 

viscosity or inflow, we have Xend = xf^m ^ ^^^n Bm,i.,min < Bm,* < Bm,i,,opt 

and Xend = 2^/^ a"^/^ when Sm,* > Bm,*,opt- However, if there is flow of wind 
material into the disk region or if there is an onset of MRI during the second 
process and viscosity becomes significant, the value of Xend will be much larger 
because of a contribution from material between X^sc and Xn^. 



2. Magnetic Field Strengths and Disk Models 

Let Bm,*,circ be the lower bound derived bv iMaheswaran Sz Cassinellil ( 19881 . 

for the meridional field to withstand the effects of meridional circulation 
near the stellar photosphere and let tdrc be the turnover time of this circulation. 
Le t Bm._*.m.tA be th e mini mum surface magnetic field required in the MTD model 
of ICassinelli et all (|2002). For Be stars with realistic rotation rates we have 
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< Bm,i,,circ < Bm,-k,mtd- Also, tfm <C tcirc, where tfm is the fill-up time 
of the pre-Keplerian disk. Table ^ gives the critical values of magnetic fields for 
the formation of disks in different stellar models. 



Table 1 . Critical values oiBm,* for MR WD and MTD type disks in different 
stellar models with specified rotation rates ijCassinelli et al2002HMaheswaranl 
|2003) 



MRWD Model MTD Model 



Spectral 
Type 


a 


(G) 


. -k , opt 

(G) 


^m.*,circ 

(G) 


^circ 

(yr) 


(G) 


03 


0.6 


437 


893 


139 


0.53 


22000 


06.5 


0.6 


101 


207 


175 


0.57 


5000 


BO 


0.6 


33 


67 


131 


0.31 


1600 


B2 


0.5 


7 


13 


25 


1.1 


335 


B5 


0.45 


1.4 


2.5 


3.6 


6.7 


64 


B9 


0.4 


0.8 


1.3 


0.5 


30.6 


30 



Here, we consider models that are appropriate for stars with different values 
of i?m,vc- Although the azimuthal component may also play a role, we do not 
include it in this discussion. 

(a) If Bm,*,min < Bm,* < Bm,*,circ, quasi-steady Keplerian disks satisfying 
the MRWD model can be formed. The surface magnetic field will be affected 
by meridional circulation near the photosphere. This can lead to changes in the 
disk over time scales of order tdrc- 

(b) When Bm,i,,circ < Bm,i, < Bm,ir,mtd, the MRWD model is appropriate. 
The surface magnetic field will be able to withstand the effects of circulation 
and persist over time periods that are long compared with tdrc- 

(c) If i?m,* > Bm,,i,,mtd, the disks will be of MTD type. If the magnetic field 
lines from the star thread the disk region and loop back to the star, like field 
lines of a dipole-type field, the angular velocity distribution in this disk will be 
the same as that in the region PinnPum of the stellar surface from which the field 
lines emerge. On the other hand, if the field lines from the star pass through the 
disk region and travel outwards to large distances from the star without looping 
back, the rotation law for th e disk will be similar to the empirical formula 
specified in the MTD model of lCassinelh et all ((2002). 



3. Discussion 



The requirement that the disk density should be positive gives the c ondition 
v'^ > {B%gj^ — B'^^^j) / S-K pyjind for disk formation ( MaheswaranI |2003 ^ . Here, 
all quantities are evaluated on the appropriate sides of a point Q on Tighock 
and Vn is the normal component of the wind velocity. Thus, for given values 
of the equatorial rotation speed and surface magnetic field strength, the wind 
velocity must be larger than a critical value for a disk to be formed. Since 
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this critical value depends on the stellar rotation speed, a faster rotation rate 
does not ensure the formation of a shock-compressed disk region unless Vn is 
correspondingly larger. 

The MRWD model does not require that the magnetic field be uniformly 
strong across the entire stellar surface. The field must have the required strength 
on the sector PinnPhm: whose spread in latitude is only a few degrees. The basic 
processes in this model do not depend on whether the field is axially symmetric. 
The results obtained for symmetrical systems will qualitatively apply to stars 
possessing magnetic fields with fi ux loops that thread the disk region. 

Observational evidence fe.g.. lDonati et alll2nnil . 1200^ : lNeineifl2nn2l . 120041 '! in- 
dicates that several stars of the typ e we consider for d isk formation are oblique 
magnetic rotators with dipole fields. iPreuss et all ()2004l l find that disk-like struc- 
tures can form along the magnetic equatorial planes of models with strong mag- 
netic fields in which the Alfven speed is faster than the rotation speed of the 
magnetosphere. Thus, when the MRWD model is applied to oblique rotators, a 
pre-Keplerian region can be formed along a plane or surface that is determined 
by gravity, centrifugal force and magnetic force acting on the disk material. 
During the second process, if the magnetic force in the disk becomes small, we 
expect that super-Keplerian material will move into Keplerian orbits in the rota- 
tional equatorial plane. In the case of oblique rotators with -Bm,* > Bfn,*,mtd, the 
results of Preuss et al imply that an MTD disk will be located in the magnetic 
equatorial plane rather than the rotational equatorial plane. 

Because the magnetic force in the disk region of the MRWD model becomes 
small at the end of the fill-up stage, th e one- armed spiral pattern of the Global 
Disk Oscillation model (e.g., Okazakii Il997l ) can be used to explain the V/R 
variability in disks of Be stars. 
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